Quantum Hall resistances of multiterminal top-gated graphene device 
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Four-terminal resistances, both longitudinal and diagonal, of a locally gated graphene device are 
measured in the quantum-Hall (QH) regime. In sharp distinction from previous two-terminal studies 
[J. R. Williams et al, Science 317, 638 (2007); B. Ozyilmaz et al, Phys. Rev. Lett. 99, 166804 
(2007)], asymmetric QH resistances are observed, which provide information on reflection as well 
as transmission of the QH edge states. Most quantized values of resistances are well analyzed by 
the assumption that all edge states are equally populated. Contrary to the expectation, however, 
a 5/2 transmission of the edge states is also found, which may be caused by incomplete mode 
mixing and/or by the presence of counter-propagating edge states. This four-terminal scheme can 
be conveniently used to study the edge-state equilibration in locally gated graphene devices as well 
as mono- and multi-layer graphene hybrid structures. 

PACS numbers: 73.43.Fj, 71.70.Di, 73.61.Wp, 73.23.-b 



I. INTRODUCTION 

Landau-level splitting in two-dimensional (2D) elec- 
tron systems under a perpendicular magnetic field reveals 
the well-known quantum-Hall (QH) effect^ When the 
Fermi energy is set between two Landau levels, a current 
circulates along the edge conduction states in a (chiral) 
direction determined by the carrier type and the direction 
of the magnetic field^^ Utilizing this chiral character 
of the edge states one can devise diverse solid-state beam 
splitters out of 2D electron gas systems£&£ 

On the other hand, due to the relativistic nature of 
the carriers, graphene, a 2D honeycomb lattice of carbon 
atoms, shows the half-integer QH effect Thus, the 
manipulation of the edge states in graphene can be of 
particular interest. However, the electrostatic deflection 
of the edge state2££i£ is not realizable in graphene due 
to the Klein tunneling of the carriers through an electro- 
static barrier. 12 Nonetheless, the possibility of controlling 
the edge-state transmission in graphene has been con- 
firmed by locally modulating the filling factor v and the 
chiral direction of the edge states i 13 i 14 Two-terminal ob- 
servation of the edge-state transmission in graphene to 
date is well explained by the complete-mode-mixing hy- 
pothesis where all edge states are equally populated at 
p-n interfaces i 13 ' 14 i 15 QH plateaus in two-terminal mea- 
surements, however, can be distorted depending on the 
sample geometry and the contact inhomogeneitiesi 16 ! 17 
Moreover, the two-terminal conductance gives the infor- 
mation only on the edge-state transmission, lacking the 
information on the reflection. Thus, for the better manip- 
ulation of the edge states of an arbitrary-shaped graphene 
device, one needs geometry-independent measurements 
that can furnish information on both the edge-state re- 
flection and transmission. 

In this paper, we report on four-terminal QH trans- 



port measurements in a top-gated bi-polar graphene de- 
vice, which show the quantization of longitudinal QH re- 
sistances as well as an asymmetry in the diagonal QH 
resistances (the meaning will be defined below). Our 
measurement scheme provides precise information on the 
reflected QH edge states in addition to the transmit- 




FIG. 1: (Color) (a) Schematic measurement configuration. 

(b) Scanning electron microscope image of the device after 
depositing a local gate, the brightest central part. The broken 
lines represent the graphene edges and the boundary of the 
PMMA insulation layer is evident by the contrast change near 
the top of the image. The white scale bar represents 2 /an. 

(c) Two-dimensional color map of Rl(Vlg,Vbg) for T at 20 
mK. A set of carrier types in the regions 1 and 2 is labeled in 
each quadrant, (d) Vbg dependence of Rl at Vlg=-2 V, 6 V 
and -7 V extracted from (c) as indicated by arrows with the 
same color (green, red, and gray, respectively). Left inset: the 
Rd (upper) and Rl (lower) as a function of Vbg at Vlg= -2 
V and 10 T, showing the half-integer QH effect. Right inset: 
a schematic top view of the device, where the local gate is 
placed on top of the colored region at the center (the region 
2). 



FIG. 2: (Color) (a) Two-dimensional plot of R l (Vlg,Vbg) 
for 10 T at 20 mK. A set of carrier types in the regions 1 
and 2 is labeled in each quadrant. (b,c) Plots of Rl(Vlg) 
extracted from (a) at Vbg=50.5 V and 22 V, as indicated 
by the white broken lines in (a). Broken lines represent the 
calculation result based on the Landauer-Biittiker formula. 



ted ones as can be obtained from two-terminal mcasure- 
mentsi 13 ' 14 Most of the results are in good quantitative 
agreement with the Landauer-Biittiker formula^* 4 - while 
unexpected resistance plateaus corresponding to the 5/2 
transmission of edge state are also observed. It may 
arise from the incomplete mode mixing and/or unusual 
QH edge states that are possibly present under the local 
gatei 18 i 19 ' 20 i 21 ' 22 i 23 ' 24 This simple four-terminal scheme 
will allow an additional insight into the edge-state equi- 
libration in bipolar graphene systems. Results of this 
study can also be utilized for the spatial manipulation of 
the Dirac fermions in graphene, which is one of the hot 
issues in graphene studies. 



II. SAMPLE PREPARATION AND 
VERIFICATION 

The sample was prepared by mechanically exfoliating 
monolayer graphene^ on a silicon substrate covered with 
a 300-nm-thick oxide layer where the silicon substrate 
was used as a back gate (BG). Electrical contacts of Cr (5 
nm) / Au (20 nm) were patterned by the method described 
elsewhere^ It was then followed by spin-coating a 30- 
nm-thick polymethyl methacrylate (PMMA, 950 K, 2% 
in anisole) dielectric layer on top of the device, which was 
cross-linked by high doses (15000 /iC/cm 2 ) of electron 
beams 2 ^ with 20 keV. Finally, a local gate (LG) of Cr 
(5 nm) / Au (40 nm) was deposited at the center of the 
device [Fig. 1(b)]. The device was cooled down to 20 mK 
in a dilution fridge (Oxford Instruments, Model AST). 
The longitudinal (Rl—Vl/I) and diagonal (Rd—Vd/I) 
resistances were measured simultaneously with two lock- 
in amplifiers, synchronized with each other at 1=2 nA 
and /=13.3 Hz [see Fig. 1(a)]. 

Figure 1(c) shows a 2D color map of Rl(Vlg,Vbg) 
measured at zero magnetic field in units of h/e 2 (all re- 
sistances will be presented in this unit afterward). In 
the figure, horizontal and diagonal bands crossing at 
(VlG)Vbg)~(-2 V,46 V) are identified, which divide the 




FIG. 3: (Color) [(a) and (b)] Two-dimensional plots of 
Rd(Vlg,Vbg) at 10 T and -10 T, respectively. A set of carrier 
types in the regions 1 and 2 is labeled in each quadrant, [(c)- 
(f)] Red curves: Rd(Vlg) extracted from (a) at Vbg=64.5 V, 
50.5 V, 35.5 V, and 22 V. Blue curves: R d {Vlg) extracted 
from (b) at Vbg=22 V, 35.5 V, 50.5 V, and 64.5 V. Thus, 
(c) and (d) correspond to the counterclockwise edge state in 
the region 1, while (e) and (f) to the clockwise edge state 
in the same region as denoted on the right. The extraction 
Vbg values are marked by the white broken lines in (a) and 
(b). Broken lines represent the calculation results from the 
Landauer-Biittiker formula. 



map into four different quadrants. These bands repre- 
sent positions of local resistance maxima arising from the 
change in the carriers in the region underneath the local 
gate (the region 2) with respect to the outside region (the 
region 1) [the right inset of Fig. 1(d)]. Figure 1(d) is the 
one-dimensional slice plot of Rl(Vbg) that is extracted 
from Fig. 1(c) at V LG =-2 V, 6 V, and -7 V. It indicates 
that the position of the secondary resistance peak, cor- 
responding to the diagonal pattern in Fig. 1(c), directly 
depends on Vlg while that of the dominant resistance 
peak, the horizontal pattern in Fig. 1(c), is almost in- 
sensitive to Vlg- This leads to a conclusion that the 
horizontal (diagonal) band is from the charge-neutrality 
point in the region 1 (in the region 2) where the carrier 
type is altered. It verifies the successful performance of 
our bipolar device i 13 ' 14 i 26 ' 27 i 28 Additionally, in the left 
inset of Fig. 1(d), we show the half-integer QH effect, 
Rd (upper) and Rl (lower), taken at Vlg=-2 V and 
S=10 T, indicating that our device consists of a mono- 
layer graphene sheet s 10 ' 11 



III. EXPERIMENTAL RESULTS 

The 2D color map of Rl{Vlg,Vbg) measured at 10 T 
is shown in Fig. 2(a). It displays several skewed blocks 
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of different resistances, implying the quantization of the 
longitudinal resistance. Details are more clearly seen in 
one-dimensional slices, Rl(Vlg), of Fig. 2(a) for fixed 
values of Vbg- hi Figs. 2(b) and 2(c), Rl(Vlg) is dis- 
played for Vbg=50.5 V and 22 V, respectively. The red 
and blue curves were taken at 10 T and -10 T, respec- 
tively. Two curves almost completely overlap with each 
other. Dominant resistance plateaus exist at zero resis- 
tance in the region -3 V<Vl G <1 V [Fig. 2(b)] and -3 
V<Vlg<0 V [Fig. 2(c)], which arise from the full trans- 
mission of the edge states when the filling factors in the 
regions 1 (z^i) and 2 (^2) are identical. In addition to 
these trivial ones, one also finds plateaus of non-zero frac- 
tional resistances such as 1/15, 1/3, and 2/3 in a certain 
range of Vlg- This quantization directly demonstrates 
that a portion of the edge states is reflected at the inter- 
faces between the regions 1 and 2 for non-identical filling 
factors v\ and 1/2, which is consistent with the previous 
two-terminal conductance measurements . 13 > 14 The quan- 
tized values of Rl are in excellent agreement with the 
calculation results following Ref. [4] as represented by 
the broken lines in Figs. 2(b) and 2(c), the details of 
which will be discussed below. 

Now, let us focus on the diagonal resistance Rd, which 
exhibits far richer features. As seen in Figs. 3(a) and 
3(b), the 2D plots of R d (Vlg,Vbg) taken at 10 T and 
-10 T, respectively, also reveal skewed blocks of differ- 
ent resistances, but with overall features much different 
from Fig. 2(a). First of all, both Figs. 3(a) and 3(b) 
show no inversion symmetry with respect to the cross- 
ing point at Vbg~46 V and V, where v\—V2=§ 
(the zero point). Nonetheless, there exists an inversion 
symmetry between Figs. 3(a) and 3(b), or equivalcntly 
a 180° rotational symmetry between them with respect 
to the zero point, which is in sharp distinction from the 
feature of Rl [Fig. 2(a)] as well as the previous two- 
terminal results i 13 ' 14 More details are revealed by the 
one-dimensional slices shown in Figs. 3(c)-3(f), which 
are again extracted from Fig. 3(a) for 10 T (red curves) 
and Fig. 3(b) for -10 T (blue curves). The correspond- 
ing values of Vbg are specified in each figure. The figures 
illustrate that the data taken at 10 T and -10 T are mir- 
ror symmetric with respect to Vlg~-1 V, which again 
confirms the inversion symmetry of Rd between the two 
opposite field directions. 

Figures 3(e) and 3(f) show the Rd measured when car- 
riers in the region 1 are holes (^i<0) for 10 T and elec- 
trons (i^i>0) for -10 T, which corresponds to the clock- 
wise edge-states in the region 1. Most of the quantized 
resistances in these figures match with the inverse of two- 
terminal conductance observed previously^ Thus, Rd in 
this region corresponds to the two-terminal (or Hall) re- 
sistance of transmitted edge states. In contrast to the 
two-terminal results, however, Figs. 3(e) and (f) show 
clear 1/2 plateaus. This indicates that the disorder ef- 
fect, which has been regarded as the cause of the ob- 
served reduction in the conductance for ^i=^2=±2 in the 
previous studyji 3 - is negligible in this four-terminal mea- 




FIG. 4: (Color) (a) Schematic top view of the edge-state con- 
figuration. Electrical contacts are shown in yellow and thick 
arrows represent the chiral direction of QH edge states. The 
case of the clockwise circulation in region 1 is shown as an 
example. We assume that K out of TV edge states are re- 
flected at the interfaces between the regions 1 and 2 (blue 
arrows). On the right, three possible equilibration processes 
in a graphene p-n-p device and corresponding values of K 
for each regime are illustrated, [(b) and (c)] Two-dimensional 
plot of the calculated Rl and Rd for a positive magnetic field 
as a function of v\ and V2- In both figures, calculated values 
of the resistance plateaus are written at the center of each 
block in units of h/e 2 . Horizontal lines represent the filling 
factor v\ where one-dimensional data in Figs. 2 and 3 are 
extracted. 



surement. On the other hand, in Figs. 3(c) and 3(d) (for 
i/i>0 at 10 T and z^<0 at -10 T), ±1/2, ±1/6, and -1/10 
plateaus are seen together with a sign change at certain 
values of Vlg- This sign change seems to be odd be- 
cause these plateaus correspond to the counterclockwise 
edge states in the region 1, where Rd is supposed to be 
negative. A careful analysis, however, indicates that the 
Rd for ui-P2<0 [the ni-p2-ni region in Fig. 3(a) and the 
Pi-112-pi region in Fig. 3(b)] is nothing but the Hall resis- 
tance for the clockwise edge states in the region 2, which 
should be positive as observed. But, still counterintuitive 
positive (1/6) plateaus appear in Fig. 3(c) even for the 
counterclockwise edge states both in the regions 1 and 2 
[^i-j/2>0 for the ni-n2-ni region in Fig. 3(a) and for the 
P1-P2-P1 region in Fig. 3(b)]. This can be accounted for 
in terms of the Landauer-Buttiker formalism for the edge 
states as will be shown below 2& 



IV. DISCUSSION 

The Landauer-Buttiker formula^ is suitable for 
studying the one-dimensional edge-state transport where 
the current through the lead a (I a ) is expressed as a lin- 
ear combination of transmission coefficients {T a p) multi- 
plied by the corresponding chemical potentials (/I/3). Ac- 
cordingly, it is required to evaluate the scattering matrix 
(S) with the elements T a ^'s and then solve the linear 
equations ( / =S~j2), which is straightforward in our case 
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because the reflection (K) happens at p-n interfaces only 
[see Fig. 4(a)]. This procedure is adopted to obtain Rl 
and Rd in Ref. [4] for the configuration that is identi- 
cal to ours. It is suggested that Rl is finite unless K 
is zero and Rd has two different values (R^ and Rd), 
depending on the measurement configurations. 
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where R a p-yS represents the resistance measured between 
the voltage leads 7 and 8 for the current injection from 
leads a to (3. In Eq. {l} it is assumed that N edge states 
circulate in the clockwise direction in the region 1, i.e., 
v\<0 (>0) for 10 T (-10 T) in our case. For a direct 
comparison with the theoretical expectation^ we adopt 
a four-terminal schematic sample configuration in Fig. 
4(a) rather than a five-terminal one corresponding to the 
actual device. Since each resistance {Rl or Rd) was ob- 
tained in a four-terminal configuration, the schematic is 
equivalent to the real device configuration, namely, Rd 
(Rl) in Fig. 1(a) corresponds to R1342 (-R4312) in Fig. 
4(a). Since -R1342 and -R4213 depend on the chiral direc- 
tion of incoming edge states Rd corresponds to R~j^ (R~^) 
for the clockwise (counterclockwise) edge states in the re- 
gion 1. Thus, Rd changes as the sign of v\ changes at a 
fixed magnetic-field direction, but it represents the same 
resistance if both v\ and the magnetic-field direction are 
altered. This explains the inversion symmetry shown in 
Fig. 3. 

Two-terminal resistance in a bipolar configuration 
studied previously-^ 3 * corresponds to i?^ only in Eq. |T]), 
which is the Hall resistance from the N — K transmitted 
edge states. It is consistent with our observation shown 
in Figs. 3(e) and 3(f), which correspond to the clockwise 
edge states in the region 1. But, the information on the 
reflected edge states is lost in the two-terminal process. 
On the other hand, R~j^ represents the eliminated Hall re- 
sistance by the K reflected edge-states, which can be ob- 
served only in a four-terminal configuration. This implies 
that the average of R^ and R D , (Rd~ R d ')/2, provides 
the Hall resistance of N incoming edge states (the sign 
for R^ is reversed because i?^ and R D are measured in 
the opposite direction with respect to the current flow). 
Thus, if more than a half of the incoming edge states are 
reflected (2K>N), R^ can change the sign, accounting 
for the odd sign of resistances seen in Figs. 3(c) and 3(d). 
On the other hand, Rl indicates the difference between 
the (two-terminal) Hall resistance of the N — K trans- 
mitted and the N incoming edge states which is always 
positive or equal to zero for K=0. In consequence, all 
the unexpected features in Rd as well as the finite value 
of Rl are well accounted for at least qualitatively. 



For a quantitative analysis, we estimated the number 
of the reflected edge states (A") in three different 
regimes^ the edge-state-transmission regime (z^r^2>0, 
N= |^i|>|^2|)j partial-equilibration regime (^i-^2>0, 
A=|^i|<|^2|), an d full-equilibration regime (v\-V2<Q, 
N=\v\ |). Based on the complete- mode- mixing hypoth- 
esis ; 13 ' 15 it is evident that K is equal to |^i|-|^2| in 
the edge-state-transmission regime, where the extra 
edge states in the region 1 (l^iH^I) are forbidden in 
the region 2. In the partial-equilibration regime, the 
excess states (I^H^il) circulate in the region 2 while 
in partial equilibration with the incoming edge states 
(\vi |). Finally, in the full-equilibration regime, the edge 
states in the regions 1 (|^i|) and 2 (li^l) circulate in 
opposite directions, propagating in the same direction 
along the p-n interfaces and equilibrate with each other. 
In the last two regimes, one can calculate K by using the 
current-conservation relationF^ 3 - Schematic edge-state 
circulation and values of the K are illustrated in the 
right panel of Fig. 4(a) for each regime. With these 
values of filling factor in Eq. (JTJ) , the Rl and Rd are 
calculated as functions of v\ and as shown below. 



(a) Edge-state transmission regime: 
Rd = 



Rl = 




"1 ^2 



(i/ x <Q at 10 T or ^>0 at -10 T) 
(i/i>Q at 10 T or z^<0 at -10 T) ' 

(2) 



(b) Partial equilibration regime: 
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Rl 
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(ui<0 at 10 T or ^i>0 at -10 T) 
(z/ x >0 at 10 T or v x <Q at -10 T) ' 

(3) 



(c) Full equilibration regime: 

f? 2| MM l1 (^i<0 at 10 T or v x >Q at -10 T) 
le^ra (^i>0 at 10 T or vi<0 at -10 T) 

h \vi\ + |f 2 | 



Rd — 



R 



(4) 



Results for a positive magnetic field (10 T) are sum- 
marized in two color maps in Figs. 4(b) and 4(c), re- 
spectively, with the calculated resistance labeled in each 
block. As seen in Fig. 4(b), unless ^1=^2, Rl becomes 
positive and is in good quantitative agreement with the 
observed results as denoted by the broken lines in Figs. 
2(b) and 2(c). The values of Rd are also well repro- 
duced by the calculation as shown in Fig. 4(c) as well 
as Fig. 3(a). Although not shown, the observed -10 T 
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data of Rd [Fig. 3(b)] are also well fit by the calculation. 
Furthermore, for v x >Q (<0) at 10 T (-10 T), the calcu- 
lated Rd becomes positive in the region where u\-V2<Q 
or |fi|>2|z/2| and Vi-f2>0, reproducing the odd sign of 
resistances found in Figs. 3(c) and 3(d). Again, all the 
observed features of Rd are well accounted for in terms 
of the Landauer-Biittiker formula. 2,4 Since all of these 
results are linked to the chiral-directional dependence of 
Rd, one can conclude that edge states in graphene also 
behave in the same way as those in a conventional 2D 
electron gas? 2 - although the QH effect itself is distinct by 
the Dirac nature of charge carriersi 9 ' 10 ' 11 

V. 5/2 TRANSMISSION OF EDGE STATES 

Most of the observed results in the color maps shown in 
Figs. 2 and 3 are quantitatively accounted for. However, 
one can distinguish some additional features which ap- 
pear as a faint diagonal bands crossing from the bottom- 
right region to the top-left region of the figures in the 
ranges Vbg>~60 V and Vbg<~30 V near the bound- 
ary where v-i changes sign. In fact, these band struc- 
tures constitute additional resistance plateaus of ~0.24, 
^0.07, and ~0.4 as indicated by arrows in Figs. 2(c), 
3(c), and 3(f), respectively. Solving Eq. Q for R L ~0.24, 
R D ~0.07, and i?£~0.4 at N=\vi\=6 gives almost iden- 
tical K values as ~3.54, ~3.52, and ~3.5, respectively. 
This indicates that these unexpected resistance plateaus 
are all produced by a single 5/2 (—N~K—6~7/2) edge- 
state transmission for |z/i|=6, |i/2|=2, and vi-V2>0, just 
before V2 changes sign or, equivalently, just before the 
edge-state transmission regime turns into the full equili- 
bration regime. This state may be related to the disor- 
der in the region 2, but no related features are present at 
2/i=i>2 = ±2 in our measurement In addition, the disor- 
der effect is not supposed to be significant in our four- 
terminal measurement. 

Another possible explanation is based on the fact 
that the plateaus take place near the boundary between 
the edge-state transmission and the full-equilibration 
regimes. One may imagine the presence of counter- 
propagating edge states: for instance, |2| + |e| clockwise 
edge states and |e| counterclockwise edge states in the 
region 2. Here, we expect that |e| should be an integer 
less than 2, because no fractional QH effect is presen t 29 ' 30 
and the maximum filling factor in the region 2 is two. 
However, a simple calculation of the total reflection K 
by adding two different K values for the clockwise edge 
states (edge-state transmission) and the counterclockwise 
ones (full equilibration) leads to |e| much larger than 2 
(about 4.7, corresponding to K=7/2), which is far from 
the expectation. 

The discrepancy can arise from the incomplete mode 
mixing for the counter-propagating edge states, which 
leads to the 5/2 transmission when only two and a half 
out of three clockwise edge states participate in the equi- 
libration process completely These counter-propagating 



edge states may originate from the compressible and in- 
compressible edge state s 18 ' 19 ! 20 formed near the bound- 
ary of the graphene by the charge accumulation as sug- 
gested for the graphene sheet placed above a global back- 
gate i^ 9 -^ However, the strange diagonal bands observed 
in our study indicate that the local gate, by some means, 
further enhances the charge accumulation. The presence 
of the QH ferromagnet state o 21 ' 22 ' 23 ' 24 can also be con- 
sidered for the counter- propagating edge states. But, 
the Zeeman energy in 10 T is only about 20 K^i which 
is much smaller than the energy difference between the 
lowest and the first exited Landau levels, ^700 K<2i Al- 
though this phenomenon can be understood by assuming 
both the incomplete mode mixing and the presence of 
counter-propagating edge states ; 18 ' 19 ' 20 ' 21 ' 22 ' 23 ' 24 further 
theoretical and experimental investigation is required for 
a conclusive account for the effect. 



VI. SUMMARY 

We studied the edge-state equilibration processes in 
a locally gated p-n-p junction of the graphene by mea- 
suring the four-terminal QH resistances. Our scheme 
enables measurements of finite longitudinal and asym- 
metric diagonal QH resistances, which furnish precise 
information on the reflected as well as the transmitted 
QH edge states. Most of our observations are quanti- 
tatively analyzed by the Landauer-Biittiker formula 2,4 
based on the complete-mode mixing hypothesis. But, 
unexpected resistance plateaus corresponding to the 5/2 
transmission of edge states are also observed. We sug- 
gest that it may arise from the incomplete mode mixing 
and/or the presence of the counter- propagating QH edge 
ctatc: . 18,19,20,21,22,23,24 

Our simple four-terminal measurement scheme re- 
ported here is not sensitive to the contact resistance and 
the sample geometry, so that it can be conveniently used 
for the edge-state manipulation of arbitrary-shaped de- 
vices such as Mach-Zehnder interferometers X Moreover, 
with the additional information on the reflection of the 
edge states, this scheme enables one to investigate the 
details of the edge-state equilibration at p-n interfaces 
such as 5/2 transmission of the edge states found in this 
study. It can be also employed to investigate the edge 
states in hybrid structures of mono-layer and multi-layer 
graphene £M1 
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